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ABSTRACT 

We report on the first flight of the balloon-borne QMAP experiment. The experiment is designed 
to make a map of the cosmic microwave background anisotropy on angular scales from 0?70 to several 
degrees. Using the map we determine the angular power spectrum of the anisotropy in multipole bands 
from £ 40 to £ ^ 140. The results are consistent with the Saskatoon (SK) measurements. The 

frequency spectral index (measured at low £) is consistent with that of CMB and inconsistent with either 
Galactic synchrotron or free-free emission. The instrument, measurement, analysis of the angular power 
spectrum, and possible systematic errors are discussed. 

Subject headings: cosmic microwave background - cosmology: observations 


1. INTRODUCTION 


ski 1992), one in each polarization. These have been fully 


Measurements of the anisotropy of the cosmic microwave 
background (CMB) are an effective probe of the state 
of the universe roughly 300,000 years after the big bang 
( White, Scott, fc Silk 1994 , Bond 1996 ). The anisotropy in 
the CMB not only provides a key test of models of struc¬ 
ture formation, but analysis of its power spectrum may 
determine values of fundamental cosmological parameters 
such as Ho, Hq, and 

The anisotropy in the CMB is mea sured with a variety 
of detector s and observing strategies (Bennett et al. 1997, 


Page 19971) . QMAP was designed to directly produce a 


degree scale map of the sky that could be analyzed for 
anisotropy. In this paper, we report on the QMAP instru¬ 
ment and the first of two balloon flights. We present maps 
of the microwave sky and determinations of the angular 
power spectrum covering scales between 9° to 0?7 in two 
frequency bands centered on 31 and 42 GHz. An analy¬ 
sis of the second flight, systematic effects and calibration 
is presented in Herbig et al. (1998, hereafter H98), and 
the data analysis methods are described in detail in de 
Oliveira-Costa et al. (1998, hereafter d098) together with 
results from combining the flights. 

2. INSTRUMENT 

The QMAP balloon-borne telescope was designed, built 
and flown twice in a period of 20 months. It consists of a 
cryogenic receiver, primary and secondary optics, pointing 
system, and data acquisition system. 

The receiver, which was cooled to 2.3 K in flight, has four 
feed horns; one at Ka-band (31 GHz) and two at Q-band 
(42 GHz) each with 6-7 GHz of bandwidth Each Ka- 


band and Q-band band feed horn has two HEMT-based 
(high electron mobility transistor) amplifiers (pospieszal- 


cha racterized and ar e described in [Wollack et al. (1997) 


and Monnelly (1996) . To account for the variable temper¬ 
ature and consequential gain drifts of the receiver, a ther¬ 
mally stabilized noise source calibrates all HEMT channels 
every 100 seconds with a 54 ms pulse. The noise power 
spectrum of each of the channels used in the final analysis 
is shown in Figure 1. The final map, discussed in Section 
4.2, is based on the total sensitivity of the instrument, 
corresponding to the inverse quadrature sum of all of the 
detector noise power spectra. 



Frequency [Hz] 

Fig. 1. — Noise power spectra of each of the channels that were 
used in the final analysis. The bottom curve is the inverse quadra¬ 
ture sum of all the channels. 

The telescope optics are similar to those use d for three 
ground ba sed observations in Saskatoon, SK (Netterfield 


et al. 1997; Wollack et al. 1997). The telescope is mounted 


on an attitude-controlled balloon platform. The beams are 
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“Though we flew a forth feed horn with a low noise 144 GHz SIS receiver (Kerr et al. 1993), it did not work on either flight. The problem 
was traced to a faulty LO in the first flight and remains undiagnosed in the ser nnri uight, iTie system n ow works in the ground-based MAT 
experiment. Also, HEMT amplifiers have improved considerably since this time (Pospieszalski et al. 1997) 
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formed by cooled corrugated feeds which under-illuminate 
an ambient temperature 0.85 m off-axis parabolic reflec¬ 
tor which in turn under-illuminates a 1.8mx 1.2m chop¬ 
ping flat mirror. The close-packed array of feeds form 
0?89±0?03, 0?66±?02 and 0?70±0?02 beams for Kal/2, 
Ql/2, and Q3/4 respectively. The beams are moved on the 
sky rapidly (~ 4.7 Hz) by the computer controlled reso¬ 
nant chopping flat which requires w 20 W to operate. The 
output of the detectors are AC coupled and sampled 160 
times during each chopper cycle. This rate is intended to 
adequately oversample the sky when the beam is moving 
at its fastest rate. The center of the 2? 7 x 2? 7 array was 
fixed at an elevation of 40? 7 for the first flight and 40? 1 for 
the second flight. The telescope is inside a large aluminum 
ground screen which is fixed with respect to the receiver 
and mirrors. 



gondola. Since the Ka/Q band receivers are relatively in¬ 
sensitive to pendulation induced atmospheric signals, no 
active elevation stabilization is required. This straightfor¬ 
ward pointing system allows us to achieve good observ¬ 
ing efficiency. During the first flight we integrated for 4.7 
hours with no interruptions or pauses. A de tailed descrip- 


tion of the pointing system can be found in Beach (1996) 


3. OBSERVATIONS 

The payload was flown on June 16, 1996 from Palestine, 
TX and November 8, 1996 from Ft. Sumner, NM. Using a 
4 mcf balloon, observations were made from an altitude of 
30 km. In both flights, the entire instrument functioned 
nominally. 

To produce high-quality maps from the data, a pixel 
should be referenced to its neighbors over multiple time 


scales, directio ns, and angular separations (Wright 1996 


Tegmark 1997). This allows one to connect each pixel to 
many others and remove inevitable instrumental drifts and 
offsets in the data analysis. Our strategy is to sweep the 
beams horizontally with a 20° peak-to-peak amplitude at 
4.7 Hz. In addition, the entire gondola is sinusoidally wob¬ 
bled, about the meridian containing the NCP, with a peak- 
to-peak amplitude of 10° with a period of 100 seconds. As 
the flight progresses, one covers the sky as shown in Fig¬ 
ure 2. In each sweep, Othrow/^beam is 20 and 30 for the Ka 
and Q-band, respectively (with Othrow/dbeam < 10, one 
has difficulty distinguishing features in the angular spec¬ 
trum). A pixel is re-observed on time scales of 0.23 s from 
the chop, 100 s from the wobble, 12 minutes between the 
upper and lower beams, and « 2 hours from the rotation 
of the earth. With a single long flight, or multiple shorter 
flights, a pixel is referenced to its neighbors in orthogonal 
directions. The chop and wobble amplitudes were reduced 
by approximately a factor of two for the second flight to 
increase the integration time per pixel. However, the basic 
scan strategy remained the same. 

4. ANALYSIS 

4.1. Data Reduction and Calibration 


Fig. 2. — Sky coverage for flight 1 (left) and flight 2 (right) with 
intensity proportional to the time spent per pixel. Breaks in the 
coverage are due to calibrations. The scan paths are shown for all 
beams. For flight 1, Q3/4 are nearer the NCP (center of plot). SK 
observed the region inside the dashed circle. 

The gondola is designed to be small, lightweight and 
rugged. This allows us to use a smaller, more manage¬ 
able balloon which increases the probability of a successful 
launch and reduces the solid angle subtended by the bal¬ 
loon. The package was recovered from both flights with 
virtually no damage and is currently being used on a sim¬ 
ilar ground-based observation in Chile 

The attitude control system obtains primary pointing 
information from either a magnetometer or a CCD based 
star camera. When in CCD mode, the centroid of the 
chosen guide star in the CCD field is found in real time. 
This information is fed into a PID loop which controls the 
torque applied to a reaction wheel at the bottom of the 


For the channels used in the analysis, there are no cuts 
to the data, no glitches were removed, and there is no de¬ 
tected interference from any on-board component ®. A to¬ 
tal of 12,884,480 temperature measurements were obtained 
in each channel. Removal of the internal calibration every 
100 seconds reduced this by 4.8%. A time-independent ra¬ 
diometric offset of 10-15 mK (peak-to-peak) was present 
in all channels during the first flight. Much of this offset 
was traced to gap in the chopper baffle. For the second 
flight the offset was 1-2 mK (part of the reduction is due 
to a factor of two smaller chop angle). In both cases the 
offset is removed from each channel as described in H98 
and d098. The absolute pointing is determined during ob¬ 
servations of Cas-A. With the pointing solution, the data 
are binned on the sky and analyzed as a map. The data 
are calibrated using Cas-A as described in H98. 


^ see tittp://dept.physics.upenn.edu/cmb.html and http://pupgg.princeton.edu/~cmb/welcome.htm] for more information about the tele¬ 
scope in its balloon and ground-based conhgurations. 

®Q1 was not used in the first flight due to interference for two hours and Q3 lost one stage of amplification. 
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4.2. Map Production 

Making a map by simply averaging the observations in 
each sky direction would be inappropriate for two reasons: 
scan-synchronous offsets would cause artifacts, and 1//- 
noise would make the map unnecessarily noisy. We there¬ 
fore adopted a more sophisticated approach, summarized 
below and described in detail in d098. 



Fig. 3.— Wiener-filtered maps from flight 1. The CMB temper¬ 
ature is shown in coordinates where the NCP is near the center of 
each map, with RA being zero at the top and increasing clockwise. 


The synchronous offsets, predominantly at the chopper 
frequency and the first overtone add mK-level signals to 
the map. This signal is eliminated by convolving the data 
with a notch filter that vanishes at these two frequencies, 
creating a data set completely independent of these offset 
harmonics and any slow variation in their amplitude. We 
choose this convolution filter to eliminate the DC compo¬ 
nent as well, and use the remaining degrees of freedom in 
the filter to render the filtered noise as close to white as 
possible. This produces a data set where each element is 
the sum of three terms: 

1. A known linear combination of the actual sky tem¬ 
peratures; 

2. A known linear combination of the higher-frequency 
synchronous offset, which we assume is constant; and 

3. filtered detector noise, whose correlations vanish for 
time separations exceeding ~ 10^ samples. 

In other words, we have a greatly overdetermined system 
of linear equations. The map we compute is the solution 
to this linear inversion problem that minimizes the map 
noise, as described in d098. For each of the five channels, 
the result of the mapmaking process is an A^-dimensional 
map vector x (containing the temperature at each of the 
N pixels) and an A^ x A^-dimensional covariance matrix S 
which characterizes the pixel noise. 


Wiener filtering is a useful tool for visualizing the data. 
It suppresses the noisiest modes in a map and shows sig¬ 
nal that is statistically significant. The Wiener-filtered 
versions of our maps are shown in Figure 3 and are given 
by Xu, = S[S -I- S]“^x, where S is the covariance matrix 
corresponding to the level of sky fluctuations that we find 
in the data. 

4.3. Map tests 

A generalized y^-test for ruling out the null hypothesis 
that a given map contains no signal, merely noise, is dis¬ 
cussed in d098. Applying this test to the Kal, Ka2, Q2, 
Q3 and Q4 channels shows that that signal is detected at 
the significance level of 17 ct, 8ct, 5ct, O.dtr and 3tT, respec¬ 
tively. To test whether this significant signal is common to 
the maps or due to systematic errors, we apply the same 
test to weighted difference maps of the form x = xi — rx 2 
for different relative weights r. The comparison of Kal 
with Ka2 is shown in Figure 4, and it is seen that they 
pass three independent tests: x is inconsistent with noise 
both for r = 0 and r = oo (in which case only Kal or only 
Ka2 are probed), but perfectly consistent with noise for 
r = 1 (which gives the difference map Kal-Ka2). Likewise, 
we find no evidence of any signal in the Q-band difference 
maps. 

To see whether the signal can be accounted for by galac¬ 
tic foregrounds, we use this test to compare the Ka-band 
maps to those in the Q-band, writing r = [r'Ka/■ 
As seen in Figure 4, this places a 2tT lower limit on the 
spectral index f3 of —1.4, which means that the signal can¬ 
not be explained by foregrounds such as free-free emission 
{(3 ^ —2.15) or synchrotron radiation (/3 ~ —2.8) alone. 
For the CMB, f3 = 0. 



0.1 1 10 -5 0 5 

Relative normalization r Spectral index P 


Fig. 4.— Evidence for common signal in the maps. The number of 
(7 at which signal is detected in the weighted difference maps xi—rx 2 
is shown (top) together with the corresponding probability that pure 
noise would give such a large x^-value (bottom). For maps at dif¬ 
ferent frequencies (right), the spectral index (3 = Inr/ \n[vKa/^Q\ is 
constrained. 


4.4. Angular Power Spectrum 



eralized eigenvalue equation Sw = ASw. When the w are 
sorted by decreasing eigenvalue A, they tend to probe from 
larger to smaller angular scales. We obtain a statistically 
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independent power estimate from each mode, then average 
these individual estimates with inverse-variance weighting 
to obtain the band power estimates shown in Figure 5. For 
more details, see d098. 


5. CONCLUSIONS 


The first QMAP flight has produced a 441 square de¬ 
gree map of the microwave sky at an angular resolution of 
~ 0?7, detecting sky signal that is statistically significant 
at > 15(7. Systematic checks confirm that the signal is 
fixed on the sky and not of instrumental origin. Spectral 
analysis shows that the dominant fluctuating component 
is the CMB. A future paper will address possible levels of 
contamination by foreground emission. T his measurement 


is con sistent with the Saskatoon results (Netterfield et al, 
19971). 
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on the left is for the Q-band and corresponds to the band power 
6Ti = l)C^/27r]i/2 = 47 + 8 ^ in a window whose mean 

and rms width is given by £ = 84 di 46. The point on the right is 
for the Ka-band and corresponds to = 49^^ pK in a window 
£ = 92 di 45. These error bars do not include a calibration error of 
12% and 11%, respectively. 
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